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Abstract. We present a collection of selected phenomena observed in Jj ¥ production 
from proton-nucleus and heavy ion collisions at energies, ranging between the SPS and 
LHC. The emphasis is placed on the related theoretical ideas or techniques, which are 
either not widely known, or offer an alternative explanation to the observed nuclear effects. 

1 Preface 

Nuclear suppression of heavy quarkonia is usually considered as a hard probe sensitive to the proper- 
' ties of the short-living medium produced in heavy ion collisions. The main challenge is to discriminate 

Q-l between initial state interactions (ISI), usually identified as cold nuclear matter effects, and final state 

q_j interaction (FSI) and attenuation of the produced quarkonium in the dense matter created in the nu- 

clear collision. While the latter is the main goal of the study, the result depends on how well can one 
single out FSI from the admixture of ISI, which cannot be measured, but only theoretically modeled. 
Important information on the cold nuclear matter effects can be learned from data on proton-nucleus 
collisions. However, there is no simple recipe for extrapolation of such information from pA to ISI in 
AA collisions. Here we identify several obstacles preventing one from doing that easily. The mecha- 
nisms of ISI considerably vary between the energies of SPS and LHC. 



2 Cold nuclear matter: pA collisions 

2.1 Evolution and absorption of a charm dipole in a medium 

Usually the nuclear ratio is evaluated with an oversimplified model [ 1| assuming that J/W attenuates 
with a constant cross section cr a b s on the way out of the nucleus, as is illustrated in Fig. Q] (left). 
Correspondingly, the nuclear modification factor has the form, 

R ^ = A^I^- e ^ TA ^' (1) 

where cr flis is treated as an unknown parameter fitted to data; T^b) — J^dzpAib^z) is the impact 
parameter dependent nuclear thickness function. The results of such an analysis at different energies 
of J I W plotted in Fig.[T](right), demonstrate a steep decrease of the effective absorption cross section 
with energy. 
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Fig. 1. Left upper: illustration for exponential attenuation of the produced charmonium. Left bottom: illustration 
for evolution of a cc propagating though a medium. Right: the break-up cross section fitted with expression (T} to 
data from fixed target experiments 1 1 1. The curves show the results of a similar fit to the theoretical calculations 
of the J I ¥ production rate with Eq. dlOt and 6=1/2, 1,2, from bottom to top respectively. 



However, as is illustrated in the left bottom part of Fig.Q] a cc dipole is produced with a starting 
small separation r Ec ~ l/m c ~ 0.1 fm, then is evolving up to the J/f mean size rj/w ~ 0.5 fm and 
eventually forms the wave function during the formation time 

2£^ / E,p 



The expansion rate can be estimated perturbatively as, 

drj Apr 

dt E e c 

Employing the uncertainty relation one gets, 



(3) 



4(0 = £ + 4 (4) 
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The dipole cross section of a cc dipole with small transverse separation rj and energy E ec can be 
approximated as ]3], 

cr- cc (E- cc ) = C(E- cc )r l T . (5) 
Then the absorption cross section as function of // W energy and path-length takes the form, 

L 

a- abs (L,E z . c ) = i j dlo- abs (F) = C(E Sc ) |^ + -^j . (6) 

o ° 

With such a varying cross section one can easily calculate the nuclear modification ratio R p a, compare 
with Eq. (Q]i, and adjust the effective cross section. The results are shown in Fig. Q] (right) by three 
curves corresponding to 6 = 2, 1, 0.5 from top to bottom respectively. Similar analysis of data |[T| 
presented in Fig.[T]agree with these calculations. 

Qualitatively, the reason for the observed falling with energy cr a b s is clear: this is a manifestation 
of the color transparency effect O. Indeed, the higher is the J/W energy, the more the initial small 
size of the cc dipole is "frozen" by Lorentz time dilation, the more transparent the nuclear medium 
is. Notice that this can also interpreted within hadronic representation as a multi-channel problem (2), 
which is equivalent, but technically more difficult description. 
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2.1.1 pA: Higher twist c-quark shadowing 



At higher energies a- a bs is affected by another time scale, the lifetime of a g — > cc fluctuation inside 
the incoming proton, which can also be interpreted as a time scale for cc pair production, 



1 



mj, x 2 m N 



(7) 



where are the usual Drell-Yan variables. This time scale is about 5 times shorter than tf. 

If t p > Ra the initial state fluctuation g — > cc leads to shadowing corrections related to a non-zero cc 
separation. This is a higher twist effect, which can be calculated within the dipole approach. In this case 
the evolution of a dipole should be treated with a more advanced theoretical tool, compared with the 
simple model used in the previous section. Here we rely on the strict quantum-mechanical description 
the dipole evolution in a medium, the path-integral technique [4|, which sums up all possible paths of 
the quarks. The evolution equation in the light-cone variables has a form of the Schrodinger equation 
for the Green function, describing propagation of a dipole from the initial longitudinal position zi and 
dipole size ri up to z 2 and r 2 B4I51 . 



d 

i— G(r 2 ,z 2 ;ri,zi) 



■77- {mi 



2 ) + V(r 2 ,z 2 ) 



G(r 2 ,z 2 ;ri,zi). 



(8) 



The real part of the light-cone potential is given by the binding cc potential, which is chosen in the 
oscillatory form, and the imaginary part is proportional to the dipole cross section, which is assumed 
to be proportional to With this conditions and realistic Woods-Saxon shape of the nuclear density, 
the calculated energy dependence of the nuclear modification factor [6 | is plotted in Fig. [2] (left). The 
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Fig. 2. Left: Nuclear ratio R pA for lead as function of cc energy, calculated with the path integral technique |6|. 
Right: Transition of the projectile gluon into a color octet dipole, g — > [cc)&, which propagates and attenuate in 
the nucleus, and interacts producing a color singlet dipole on a nucleon with coordinates z, b. 



primary rise of R p a occurs in the same energy range, which was presented in Fig. Q] This rise is due 
to the color transparency effect, and it would continue as is depicted by the dashed curve, if no other 
effects were contributing. However, the charm shadowing effect, illustrated in Fig. [2] (right), which 
onsets at long production time Eq. causes an additional sizable suppression, and the total result is 
plotted by solid curve. 
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2.1.2 pA: Charmonium suppression at RHIC/LHC 

At the energies of RHIC and LHC all coherence time scales become long, tj > t p » R A , so the Green 
function Eq. © approaches the asymptotic limit of a "frozen" dipole size, G(r2, Zi\ t\, Zi) = <5(i"i - Vj) 
for zi - Z\ ~ Ra- In this case the path integral formalism is essentially simplified, and the nuclear 
modification factor at a given impact parameter and longitudinal coordinate z of the color exchange 
interaction, takes the form J6) , 

R P A(b, z) = J d 2 r T K (m c r T ) r\ ^( rr)e 4^)r-(M4^)r + (M (9) 

where T-{b,z) = £ m dz'p A (b,z'); T + (b,z) = T A (b) - T-(b,z); T A (b) = T-(b, oo); and 

9 1 

o-ccg(r T ) = -o- ec (r T /2) - -<T Ec (r T ) (10) 
4 8 

Naively, one could think that a color octet cc pair propagating through the nucleus (see Fig. [2] right) 
can experiences color exchanges remaining in the color octet states. So it cannot be absorbed, i.e. does 
not attenuate, and should not lead to initial-state shadowing. However, every process with a nonzero 
cross section shadows itself with this cross section. The cross section of cc production is given by the 
three -body dipole cross section o- ecg (7), this is why it enters the exponent in Eq. (O |6|. Both cross 
sections cr- CCCJ and cr ec steeply rise with rapidity <x S£ . cc Q^ixj) cc e °- 2&& i^ according to the parametrization 
El fitted to DIS data from HERA. 

The results of calculations are compared with RHIC data (9) at Vs = 200 GeV in the left panel 
of Fig. [3] The parameter free calculation IflOl well agree with the data. The gluon shadowing correc- 
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Fig. 3. Left: Data |9| for the p T -integrated nuclear suppression factor RjAuiy) for J I W produced in d-Au collisions 
with rapidity yaX^[~s = 200 GeV. The upper solid curve presents the results of calculations 1 10] at ->fs = 200 GeV. 
Gluon shadowing is negligibly weak I11I5I12I . The lower solid (dashed) curve shows predictions for proton-lead 
collisions at -\fs = 5.5 TeV, including (excluding) gluon shadowing. Right: All curves have the same meaning as 
in the left panel, but for T production. 



tion was found to be negligibly small within the shown range of rapidity, because the gluon radiation 
coherence length is too short. This is confirmed by the global analysis Ifl2l and dipole model calcula- 
tions [5]. Notice that these data are frequently described within a kind of an "upside-down" scenario. 
Namely, the break-up dipole cross section, which depends on energy and is well known from HERA 
data, is assumed to be an unknown constant. However, the magnitude of gluon shadowing, which is 
currently a controversial and model dependent issue fl3l . is assumed to be known. 
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Predictions for J/¥ suppression at the LHC energy V* = 5.5 TeV is also shown in the same plot. 
The dashed curve includes initial-state c-quark shadowing and final-state dipole break-up, given by 
Eq. (0. Solid curve differs by inclusion of gluon shadowing, which is visible, but still weak. 

Analogous mechanisms of nuclear suppression of heavy quarkonium production are applicable 
also to T. Suppression of the radial excitations T(2S), T(3S) is expected to be similar to T (compare 
with ¥' vs J I ¥), since it is mainly controlled by the size of the produced heavy dipole. Our predictions 
for the energies of RIHC and LHC are plotted in the right panel of Fig. [3] 



2.1 .3 pA: why nuclear suppression scales in xf 

It has been observed in fixed target experiments 0141151 that nuclear suppression of J/ ¥ produced 
at forward rapidities in pA collisions are strongly suppressed. Moreover, the Xf dependence of the 
nuclear suppression factor is nearly independent of energy. Since large xf correspond to small X2 of 
the target gluons, it is tempting to relate the increase of nuclear suppression with the coherence effects, 
like gluon shadowing 1 6 1 . In this case, however, one would expect the suppression factor to scale in Xi, 
rather than Xf. Remarkably, the an enhanced suppression of particles produced at forward rapidities 
has been observed in any process, hard or soft, studied experimentally. A natural explanation was 
proposed in lfl6l . 

Multiple interactions of the projectile hadron and its debris propagating through the nucleus should 
cause a dissipation of energy, what should result in a deficit of energy in a nuclear process at large 
Xp. This intuitive expectation is supported by consideration of the Fock state decomposition. The 
projectile hadron can be expanded over different states which are the fluctuations of this hadron. In 
the limit of infinite momentum frame those fluctuations live forever. One can probe the Fock state 
expansion by interaction with a target. The interaction modifies the weights of the Fock states, a nuclear 
target enhances higher Fock components in the projectile hadron. In each Fock component the hadron 
momentum is shared by the constituents, and the momentum distribution depends on their multiplicity: 
the more constituents are involved, the smaller is the mean energy per a constituent parton, i.e. the 
softer is the fractional energy distribution of a leading parton. So on a nuclear target the projectile 
parton distribution falls at xf — > 1 steeper than on a proton. Apparently, this effect scales in xp- 
Further details of this mechanism of suppression and numerical results can be found in 11 161 1 31 . 



2.2 pA: Cronin effect 

Nuclear targets modify the transverse momentum distribution of produced particles, suppress it at 
small but enhance at medium large pj. This effect named after Cronin, can be calculated within the 
dipole model 111711 8 1 . A simple description of the /^-dependent nuclear modification factor for J/¥ 
production was proposed in lfl9l . Available data on J/ ¥ produced in pp collisions are well described 
by the following parametrization of the p j -dependence {pj < 5 GeV), 



dcr pp (J/¥) 
dph 



1 + 



p\ ' 
^pl), 



(ID 



Following 1191201 . the simple way to calculate the nuclear modification factor is to make a shift in the 
mean square of Jj ¥ transverse momentum on the nuclear target {p T ) pA = {p\)pp + A pA . 



R pA (p T ,b) = 



(pl)R 



pA 



1 + 



P\_ 

6{P T ), 



1 + 



6[{p 2 T )+A pA (b)]) 



(12) 



{p 2 T )+A pA (b) 

Here A pA {b) = (p T ) pA - (Pj)pp is the broadening for a J/¥ produced at impact parameter b. It was 
calculated in CD within the dipole approach as, 



= o %°-dip{r T )\ _ A T A Q>\ 



r T =0 



(13) 
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in good agreement with available data for J/ ¥ broadening 0221231 . 

Formula ( TT2b is compared with data Il24l in the left pane of Fig. [4] and predictions for RHIC and 
LHC are presented in the right pane. For the energy dependence of we use the parametrization from 
(251, (p 2 T ) = [-2.4 + 0.6 In s] GeV 2 . 




P T (GeV) p T (GeV) 

Fig. 4. Left: The exponent a = 1 + ln(R pA )/ ln(A) as function of p T calculated with Eq. < 12b in comparison with 
data from the E866 experiment [24]. Right: Prediction for the p^-distributions of J/W produced with rapidities 
y = 0, 2 in p-Au collisions at -y^ = 200 GeV and in p-Pb collisions at = 5.5 TeV. 



3 Initial state interactions (ISI): transition from pA to AA 

The nuclear modification of J I ¥ production in AA collisions originate from the initial (ISI) and final 
state interaction (FSI) stages. The former includes interactions during the propagation of the nuclei 
through each other, while the latter corresponds to the interaction with the created matter, which occurs 
at a much longer time scale, when the high-energy parts of the nuclear debris are already far apart. 
Such a factorization into the two stages is possible, because no interference between them is possible. 

We consider here several important aspects of ISI. Although it is tempting to extrapolate our ex- 
perience with pA interactions to ISI in AA collisions, such a procedure is not as straightforward as it 
looks at the first glance ifTTl , In particular the so called "cold nuclear matter" in AA collisions, turns 
out to be not cold at all. 



3.1 Broadening of ./ /'/' in pA and AA collisions 

Broadening is predominantly an ISI effect, since it is not affected by FSfl. This is a sensitive and un- 
biased probe for the properties on the "cold nuclear matter". Indeed according to Eq. ( TT3l broadening 
is proportional to the medium density integrated along the parton path length. Thus, one should expect 
an universal broadening effect in pA and AA collisions, provided that the total path lengths are equal. 
However the comparison done recently in [26] and presented in Fig.[5]does not confirm such an uni- 
versality. Data show that broadening in nuclear collisions is about twice as large as in pA. This is a 

1 The cross section of J/¥ absorption in a medium is expected to rise with energy, i.e. with p T , what should 
lead to shrinkage of the p T distribution. We found this correction rather small. 
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Fig. 5. Left: The mean transverse momentum squared of J/ ¥ produced in pA and AA collisions as function of the 
mean path length in the nuclear medium. Data points are from fixed target experiments 1261 at 158 GeV. Right: 
time - longitudinal coordinate plot for charmonium production in a nuclear AB collision in the cm. The solid and 
dashed lines correspond to the nucleons and radiated gluons trajectories, respectively. 



clear evidence of a significant difference between what is called "cold nuclear matter" in AA collisions 
and that in pA. 

In fact, such a difference was predicted in [27 28 1, where it was noticed that both the projectile 
gluon and the produced J/ ¥ (or the cc dipole) interact with bound "cold" nucleons in the case of a 
pA collision, as is illustrated in the left pane of Fig.Q] However, in AA collisions those target nucleons 
have already had a chance to interact with other nucleons in the beam nucleus, so they are the debris of 
the colliding nucleons and are not "cold" any more. Apparently J/¥ should interact with this nucleon 
debris with an increased cross section. In particular, such excited nucleons should be accompanied by 
radiated (on-mass-shell) gluons, as is illustrated in the right pane of Fig. [5] 

The radiated gluons participate in the J I ¥ break-up, as well as in broadening. Each on-mass-shell 
gluon contributes as cr / - f 0\u* ^o~™. The mean number of radiated gluon per nucleon, (n g ), can 
be estimated relying on the time scale for gluon radiation, P, = 2 E q x(l - x)/(x 2 m 2 + k 2 ), where x is 
the fractional momentum of the radiated gluon. 



We see that at the energy of SPS every participating nucleon has about one extra gluon, increasing 
its interaction cross section. However, the amount of such gluons steeply decreases with the collision 
energy and vanishes at the energies of RHIC and LHC. 

The results of calculations B27I281 are compared Fig. |6]with data on J/ ¥ suppression in minimum 
biased events (left) and as function of centrality (right). Apparently, the results with no radiated gluons 
grossly overestimate data for heavy nuclei and for central collisions, what lead to a conclusion about 
anomalous J/¥ suppression. At the same time, the data are well explained with (n g ) ~ 1. These 
result point at a possibility that the so called "anomalous" suppression of Jj ¥ observed in heavy ion 
collisions at the SPS, does not signal about a hot matter produced in final state, but is mainly a result 
of enhanced ISI. The used as a baseline simple model assuming that ISI interactions of J/ ¥ are the 
same as in pA collisions, is just incorrect. This also explains why no jet quenching has been observed 
at SPS 021 . 

Increase of the dipole-Af* cross section by factor 1 + 0.75(n g ) * 1.75 also well agrees with the 
observed larger broadening of J/¥ in AA compared with pA collisions, as is depicted in Fig. [5] (left). 
Since FSI does not affect broadening, the observed strong increase confirms that both "anomalous" 
effects of enhanced J/ ¥ suppression and broadening originate from ISI. 




(V^ = 20GeV) 
(Vs = 200GeV) (14) 
(Vs= 1200 GeV) 
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Fig. 6. Left: Ratio of to Drell-Yan cross section as a function of centrality in S-U collisions at 200 GeV/c, with 
(n g ) = 0, 0.5, 0.75. The data points are from 1301 . The calculated curves are normalized at E T = to the ratio 
observed for p-p collisions. Right: The J/ if/ total cross section as a function of the product AB of the projectile 
and target atomic mass numbers at 200 GeV/c, for (n g ) = 0, 0.5, 0.75. The data are from 1311 



3.2 Mutual boosting of the saturation scales in AA collisions 

The partonic structure of a hadron is known to depend on the hardness of the probe, the higher is the 
resolution, the more partons at small x (and the less partons at x — > 1) is resolved. This is controlled by 
the DGLAP evolution. A nuclear target provides a harder probe for the partonic structure of the beam 
hadron compared with a proton target, because partons propagating through the nuclear target get an 
additional transverse kick, known as broadening. Therefore, the projectile proton in pA collisions ac- 
quires more partons at small x than in pp interactions. In the case of AA collisions all the participating 
nucleons in both nuclei change their partonic structure acquiring more low-x gluons. This leads to an 
increase of the interaction cross section for such nucleons, resulting in enhanced broadening, which in 
turn excites the partonic structure of nucleons from another side even more. Such a mutual boosting 
of parton density in the colliding nuclei leads to a rise of the saturation scales in the colliding nuclei, 
compared with pA collisions. This process is described by the following bootstrap equations [33 1 

3tt 2 

Q]b(Xb) = —a s {Q] A + Ql)x B gN{x B , Q] A + Go) T B (15) 
3tt 2 

Q 2 sA (x a ) = —a s {Q] B + Ql)x A g N (x A , Q] B + Q 2 ) T A . 

Here Q 2 A (x A ) is the boosted saturation scale in the nucleus A; Q? is an infrared cutoff providing the 
correct behavior in the soft limit (see in 13310 : T A ^ are the nuclear thickness functions of the colliding 
nuclei. 

Solving these equations one can calculate the ISI suppression of Jj ¥ in AA collisions and compare 
with the results obtained with "normal" value of Q s , the same as in pA collisions. This comparison 
presented in Fig. [7] (left). While at the energy of RHIC the boosting effect is rather mild, in the LHC 
energies the boosted saturation scale makes the nuclear medium significantly more opaque, and J/ W 
is several times more suppressed compared to the simplified expectations. 

The rise of the saturation scale means an increase of broadening Il22l in AA versus pA collisions. 
This effect looks similar to what was observed at SPS as we discussed in the previous Sect. 13.11 
However that one disappears at the energies of RHIC and LHC, while the boosted saturation scale is a 
coherence effect, which sets on only at high energies. 

Thus, we again conclude that the so called "cold nuclear matter" in AA collisions is not cold. 
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Fig. 7. Left: The nuclear ratio R AA for central (b = 0) collisions of gold-gold (RHIC) and lead-lead (LHC) vs 
impact parameter r. The upper and bottom dashed curves correspond to calculations with the same saturation scale 
as in pA at *>fs = 200 GeV, 5.5 TeV respectively. Solid curves are calculated with the boosted saturation scale, 
which makes the nuclei more opaque. Right: The T-integrated RAA(b) for J/W suppression by the ISI. Solid and 
dashed curves present the results at y = including and excluding the effect of double-color-filtering, respectively. 
The upper and bottom pairs of curves (solid and dashed) correspond to y = and energies -\fs = 200 GeV and 
5.5 TeV respectively. The dotted curve is calculated at y = 2. 



3.3 Double color filtering 

In a Glauber-like approach the survival probability of a cc dipole propagating through the colliding 
nuclei at the ISI stage is a product of the survival probabilities in each of the nuclei. This also seems 
to go along with the conventional intuition. We demonstrate, however, that at high energies this is not 
correct, and the double-color-filtering effect makes the transition from pA to AA not so straightforward, 
as usually believed IfTTI . 

As an example, let us take a piece of nuclear matter of constant density p. Classically, the survival 
probability of a dipole propagating a path length L in the medium is P(L) = e~°~" bspL (compare with 
(Q])). However, the absorption cross section is an average over the dipole size weighted with the dipole 
wave function, cr a i, s = (ctqq^t)) = J <PrTWqq(TT)?'o~qq(TT)- At high energies If » L the whole 
exponential must be averaged, rather than just the exponent (like in the Glauber model), 



P(L) = ({^H = — 

X ' 1 + (Tabs 



pL 



(16) 



For the sake of simplicity we assumed a gaussian form for W qq (rT), and (Jqq(rj) — C r*. The result 
shows that the effect of color transparency makes the medium more transparent (as expected). 

Naively, one could guess that the survival probability for simultaneous propagation through the 
two nuclei with the path lengths La and Lb has the form of a simple product, 

P(L A , L B ) = P(L A )P{L B ) = — —J— — . (17) 

(1 + Cr abs pL A ){ \ + (TabsPLB) 

However, in this case one should average over rj the product of two exponentials with the result, 

P(La,L b ) = — 1 (18) 

1 + o- abs p(L A + L B ) 

which is quite different from (TTTI i. This difference comes from the effect of double color filtering, 
which makes the medium more transparent because filtering in one nucleus reduces the average size 
of the survived dipoles, so that the other nucleus becomes more transparent. 
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Numerically, this effect is not very strong, as is demonstrated in Fig. [7] (right). It works in the op- 
posite direction to the boosting effect, which is numerically stronger and makes the colliding nuclei 
more opaque. The combined effect of both phenomena on the J/ W suppression at the ISI stage (in- 
cluding also the Cronin effect) is shown by solid curves in Fig.[8](left), We see that the boosting effect 




P T (GeV) p T (GeV) 

Fig. 8. Left: Nuclear ratio for J I W produced in central (b = 0) AA collisions including only the ISI effects. Dashed 
lines correspond to (from top to bottom) Cu-Cu and Au-Au at *ifs = 200 GeV, Pb-Pb at VI = 2.76 and 5.5 TeV, 
respectively. The solid curves include the effects of boosted saturation scale and double color filtering. Right: 
RHIC data I34I35I36I for J/W production in central collisions of Cu-Cu (closed circles and squares) and Au-Au 
(open circles) at ^ = 200 GeV and ;/ = 0. The curves are the same as in the left pane of this figure, but corrected 
for the FSI effects, calculated with the transport coefficient parameter for gold-gold q = 0.6 GeV 2 / fm adjusted 
to the data in |20|. Raa at the LHC energies is predicted with go = 0.8 GeV 2 / fm, extracted in 1 38 ] from data for 
nuclear quenching of high-p r hadrons [39| (see text). 



significantly enhances the J/ W production rate at large pj at LHC. This happens because broadening 
is boosted. 



4 AA collisions: combining the ISI and FSI effects 

Considering Jj ¥ produced with where Data from RHIC for the pj dependence of Jj W production 
with a reasonable accuracy are currently available at pj < 5 GeV. In this kinematical domain one can 
evaluate the characteristic time scales. The production and formation times in the rest frame of the 
medium are given by, 

t' p = 1 < 0.07 fm (19) 



f f = — < 0.5 fm. (20) 

are shorter than the time scale of medium creation, to ~ 1 fm, and much shorter than the mean path 
length in the medium L ~ 5 fm. So we conclude that differently from the high energy limit, where a 
"frozen" cc dipole propagates through the final state medium, in this case it is a fully formed Jj W, so 
one can rely on the Glauber approximation. 

The dipole cross section is related by Eq. (TT3T > to parton broadening, i.e. to the transport coefficient 
q, which is defined as the broadening rate per unit of length, q — 2p dcr(rT)/drj\ r= o. Thus, the dipole 
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break-up rate is also related to q as dS (rj, l)/dl = -^qr\ \ 191201 . The FSI modification factor gets 
the form, 



R(s,p T )= - j d(f> exp-i<r 2 > jdlq(s + l) 



(21) 



For the coordinate and time dependence of the transport coefficient one can employ the popular 
model Ifl9l20l 

„ A qo to n part (b, t) 

q(b, r, t) = — — — (22) 

t n par ,(0, 0) 

The magnitude of the transport coefficient was adjusted in II 191201 at q = 0.6 GeV 2 / fm to repro- 
duce the data, as is depicted in Fig. [8] (right). The ISI effects are included. The smallness of qo resolves 
the puzzle of the unexpectedly large (20 times larger) value extracted from jet quenching data within 
the energy loss scenrio 071 . 

The nuclear modification factor at the energies of LHC was also predicted in ifTOl relying on 
the result qo = 0.8 GeV 2 / fm of the analysis [38 1 of ALICE data for high-p^ hadron suppression 
1 39 1 performed within a scenario based on dipole attenuation . Notice than the recent more advanced 
analysis ll40l of RHIC and LHC data for high-p^ hadrons employing the path-integral technique, led 
to close, but somewhat higher values qo - 1 .6 and 2 GeV 2 / fm respectively. The difference is probably 
related to the approximations made above. In particular, the quadratic rise of the dipole cross section, 
cr(rr) oc rf up to a size of J/ ¥ certainly overestimates the absorption cross section in Eq. (fJTJ. Also the 
formation time Eq. (l20b was completely neglected, but it would have been more accurate to use the path 
integral technique. Both this approximations lead to to a more opaque medium, i.e. to a diminished 
value of qo adjusted to data. A more comprehensive analysis of 7/ ¥ suppression is in progress. 

Another potential source of a missed dynamics, which may lead to a reduction of the effective qo is 
the thermal mechanism of J/ ¥ production |41 42 43 1 due to coalescence of charm quarks abundantly 
produced at high energies. Although, no clear signal of this mechanism has been ob served in the 
experiments at the SPS and RHIC, it seems to be the only way to understand the unusual behavior 
of the nuclear effects as function of pr and centrality observed at the LHC 11441 . A more detailed 
discussion of this mechanism is beyond the scopes of this paper. 



5 Summary 

- Interplay between the effects of color transparency in the final state and charm shadowing in the 
initial state leads to a peculiar energy dependence of the effective absorption cross section of J/ ¥ 
produced in pA collisions. It rises up to some energy, as is confirmed by fixed target experiments, 
then it starts falling with energy being suppressed by charm shadowing. 

- At the energies of RHIC and LHC both effects reach the asymptotic regime and become the main 
source of J/¥ suppression, although they are higher twists. The leading twist gluon shadowing is 
found rather weak. 

- The transition from pA to AA collisions is not as straightforward as is usually believed. The pro- 
duced cc dipole interacts not with bound nucleons, but with the colored debris of those who have 
already had a chance to interact prior meeting the charm dipole. This leads to an increase break-up 
cross section of the dipole, as well to an enhanced broadening of the primordial projectile gluon. 
Both effects, anomalous (compared with pA) J/ ¥ suppression and broadening were observed in 
AA collisions at SPS, and both vanish at higher energies of RHIC and LHC. 

- New effect of a boosted saturation scales affecting ISI in colliding nuclei onsets at high energies. 
Although this effect has a different origin, it acts similar to what was observed at the energies of 
SPS, namely the increase of the saturation scale leads to a stronger suppression of J/¥ by ISI 
and to a larger broadening. While the magnitude of the effect is rather mild at RHIC, it is grossly 
enhances at the energies of LHC. 

- The affect of double color filtering acts in the opposite direction to the boosting, it makes the nuclei 
more transparent. 
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-J/W production offers an alternative probe for the transport coefficient of the medium created 
in heavy ion collisions. The final state attenuation of J/W is controlled by the same transport 
coefficient as parton broadening and energy loss. The small transport coefficient found from the 
analysis of Jj W data is close to what has been predicted, but is much smaller than the result of jet 
quenching analyses based on the energy loss scenario. 
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